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Abstract. The phenomenon of the so called Fermion condensation, a phase 
transition analogous to Bose condensation but for Fermions, postulated in the past to 
occur in systems with strong momentum dependent forces, is reanalysed in a model 
with infinite range interactions. The strongly non Fermi Liquid behavior of this 
system is demonstrated analytically at T = and at T ^ in the superconducting 
and normal phases. The validity of the quasiparticle picture is investigated and 
seems to hold true for temperatures less than the characteristic temperature Tf of 
the Fermion condensation. 



1 Introduction 

Recently, proceeding from the Landau approach [1] which treats the ground state 
energy E of a Fermi system as a functional of its quasiparticle distribution n P)CT , 
Khodel and Shaginyan [2] proposed the idea of the so called Fermion Condensation 
(FC). This phenomenon can occur in systems with a sufficiently strong momentum 
dependence of the effective interaction between particles giving rise to a downswing 
of the quasiparticle dispersion in the vicinity of the chemical potential /i (Fig.l). 
A distinct feature of this phenomenon is the emergence of a flat (dispersionless) 
portion of extension f2 in the quasiparticle spectrum £ P)CT = ^p-, located at the 
chemical potential fi: £ P)CT = //, for p e fl This plateau is obtained from the 
Maxwell construction as shown in Fig.l and as explained in [2]. As a consequence of 
the high degeneracy of this spectrum the distribution n p (T varies with p continuously 
with a finite derivative dn p>a /dp in contrast to np{p) = 6{p — pp) of ordinary 
Fermi liquids. A macroscopic part of quasiparticles can then occupy the region fl 
( the Fermion condensate), a feature very similar to Bose-Einstein condensation of 
bosons. The condensate wave function turns out to be degenerate. This degeneracy 
may be removed by the interactions excluded from the initial functional E[n P}(7 ], e.g. 
pairing correlations [2]. As a result, the plateau in the spectrum of single-particle 
excitations is distorted, but, as a rule, the spectrum remains quite flat since the 
strength of pairing force is rather small. 

Some difficulties with the concept of FC have later been discussed by P.Nozieres 
within the Hartree-Fock approach by using a schematic model of infinite range forces 
[3] providing an extremely strong momentum dependence of the effective interaction 
entering the Fock term. In fact, the radius r s of the force is supposed to greatly 
exceed the interparticle distance r and then the model [3] provides us with the 
leading terms of the r /r s expansion. This model proves to be very useful, since it 
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enables us to analytically evaluate most of the basic quantities. We will adopt it to 
a large extent in this paper. 

The purpose of the present work is to further elucidate the rather unusual prop- 
erties of Fermi systems which undergo FC. The outline of the paper is as follows. 
In Sect. 2 we will address the question of the interplay between superconductivity 
and FC at T = 0. Sect. 3 is devoted to finite temperatures, and in Sect. 4 we discuss 
the validity of the quasiparticle picture. Finally in Sect. 5 we conclude our paper. 



2 Interplay of superconductivity and Fermion con- 
densation. 

In order to formulate our studies in a convenient way we take up the model of infinite 
range forces [3] in a slightly more general way. We start out from the following 
Hamiltonian 

H = H + Hi nt , (1) 

where 

Hq = £ Pi(7 £ p ap )(T ap i(J , £ p = £ p — Hi (2) 

and 

Hint — 2^Pi,P2,q,aiv,04[^o( ( l)^ai,a3^Q2,Q4 + ^' ^s( ( l) CT ^a 1 ,a 3 "a2,O4] a pi,Qi a P2,a2 a Pi+ c l' a 3 a P2-q,a4 

(3) 

with e p = p 2 /2m, U (q) = U 5(q) and U s (q) = U s S(q). The constants U and U s 
are assumed to be small compared to the Fermi energy e° F = p 2 F /2m. The only 
difference with the model [3] is that the Hamiltonian (1-3) contains an additional 
spin exchange mixture. 

In order to study the superconducting properties of our model let us calculate 
the ground state energy from the usual singlet BCS wave function 

\BCS >= Tl(u p + Vpd+^atp^ \vac > . (4) 

This ground state is known to be the vacuum to the quasiparticle operators a P)(J (i.e. 
a p ^\BCS >= 0) which are related to the original operators through the Bogoliubov 
transformation 

o P ,+ = u P a p,+ + v p a± p _, a p _ = u p a p _ - v p a± p + , (5) 
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with the coefficients u v and v p obeying the normalization condition u 2 +v 2 = 1. The 
expectation value of H is easily calculated and we obtain 

F =< BCS\H\BCS >=E-nN = E p [2^n p + V x n\ - V 2 K* p K p ], (6) 

where 

Tip = 2^°" *^ a p,o- a p,CT ^ = k p = ^ a p,+ a -p,— > = UpVp, (7) 

are the occupation numbers and the pair distribution function, respectively. A 
straightforward calculation yields V\ = —(U + 3U S ), V 2 = 3U S — U . For the model 
[3] with U s — 0, the condition V± > 0, necessary for the existence of the Fermion 
condensate, is fulfilled only for attractive forces Uq < that also leads inevitably to 
Cooper pairing, since in this case V 2 > 0. However, by adding a spin-dependent term 
U s < we can, for example, take repulsive Uq > and choose the spin-dependent 
constant U s < such that FC exists while Cooper pairing not. Note, that the 
Hartree term reducing to a renormalization of the chemical potential \i is omitted in 
(6). We also want to indicate that the infinite range of V\ and V 2 only is a convenient 
tool to trigger FC and to qualitatively investigate the interplay between pairing and 
FC. 

At T = 0, the BCS extremum equation = 0, corresponding to the minimum 
of E provided V\ + V 2 > 0, has the well known form(see, e.g. [4]): 

Z P -A p (u 2 p -v 2 p )/2u p Vp = 0, (8) 

where £ p = e p — /i and the Landau quasiparticle energy e p = J^- and the gap 

a p = are s iven b y 

e p = e° p + V inp , (9) 

and 

A p = V 2 u p Vp. (10) 

It is worth noting that the substitution of the Fermi step n p = 8(p — pp) into eq.(9) 
yields a downswing in e p with a vertical slope lying exactly at the Fermi surface. 

Upon substitution of (10) into (8) with the replacement of u 2 — v p by 1 — 2n p we 
obtain 

£ P - V 2 (l - 2n p )/2 = 0. (11) 

Together with (9) this yields 

H + V 2 /2-p 2 /2m p 2 /2m- fi + Vi/2 



n„ = 



V 1 + V 2 ^ V 1 + V 2 



e P = r' li ~ v 2 . (12) 
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From (12) we see, as discussed in the introduction, that e p — /i is dispersionless if 
V2 is zero, i.e. for vanishing pairing interaction. 

The occupation numbers n p must lie between and 1. We therefore have n p = 1 
for p < Pi and n p = for p > pj. Upon inserting these conditions into (12) one finds 

p) = 2m(n + V 2 /2), p\ = 2m{n - V 2 /2 - Vi). (13) 

With the help of (13) n p and k p can be rewritten in a very compact form 



n 



P}-P 2 ^(P}-P 2 )(P 2 -P1) 

V-2 2> K P=~ 2 2 ' Pi<P<Pfi ( 14 ) 

while outside the above domain n p coincides with the usual Fermi gas distribution 
(see Fig. 2) where k p is zero. We infer that the minimum of the free energy F is 
attained inside the functional space [n p ] only if the momentum p belongs to the 
condensate region, otherwise this minimum gets into a boundary point of the space 
[n p ] at which the gap A p vanishes. 

An additional relation, linking pi to pf and appropriate for the evaluation of 
the chemical potential /z, follows from the equality between the quasiparticle and 
particle numbers. We have 

p 3 F _ pf [Pf dpp 2 pj - p 2 



f'f>f dpi 

J Pl 7T 2 



3tt 2 3tt 2 ' J n vr 2 p)-pf 
The result of the integration is 

2 (j}y P^) 

Hp} -pf) 

In this article we shall concentrate on the case V2 << V\ in which the pairing forces 
are weaker than the particle-hole interaction V\ since the case V\ ~ V2 has already 
been analyzed in [3] while the choice V\ « V 2 leads to ordinary BCS theory. The 
quasiparticle spectrum E p can be evaluated by rewriting (11) as 2n p = 1 — 2£ P /V2 
and comparing to the standard BCS form 

2n p =\-i P jE p (12'), 

where E p is the quasiparticle energy. Then we find that the spectrum E p in the FC 
region is dispersionless as in [3]: 

e p = v / € + A P = y' Pi<p<Pf- (16) 
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Relation (16) can also be verified directly using (7,10,12,14). With (16) it is straigh- 
forward to show that the usual expression for the pair distribution 

is compatible with (14). 

In the case V 2 = 0, the quasiparticle distribution n p (12) transforms to 



n 



ji — p 2 /2m 



P-<P< P° f , (17) 



where p® = ^2m/i and p° = ^2m(/i — Vi).In this case the expression (17) coincides 
with that determined by the minimum equation [2]: j^- = fj, and the spectrum e v 

from (12) does have the plateau e p — n inside the interval p® < p < p®. 

Using (13), (15) and keeping only leading order terms in Vi/e° F we get for the 
correction to the chemical potential: 

pL - pLF ^-(V 1 + V 2 ) 2 /A8pi F ). (18) 

Here we have used the formula fXp — p F /2m + Vi/2 for the HF chemical potential. 
With the chemical potential (18) one can evaluate the spectrum of the collective 
excitations that are sound vibrations. The sound velocity c s is found from the 
relation mc 2 s = and is given by 

° s 3mA + 48(/i F ) 2 )' 1 j 

We see that the rearrangement of the single particle degrees of freedom influences 
little the collective ones since the second term is a small correction to the first one. 

From the above considerations we can draw our first conclusions. For example, 
let us consider the limit V 2 — > 0. From (12, 12', 16) we see that the distribution n p 
is an absolutely smooth and continuous function of V 2 down to V 2 — 0. The fall 
off width of n p is practically independent of V 2 and is dominated by the particle- 
hole component V]_. Furthermore, the superconductivity order parameter (14,14') 
is also a continuous function of V 2 down to the limit V 2 — which thus differs 
from zero even if A p — 0. In other words, we can say that even in the absence 
of pairing forces, if the situation for the Fermion condensation prevails, the system 
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spontaneously undergoes a transition to the state with broken gauge invariance that 
is traditionally associated with superconductivity. This is also clear from the fact 
that, according to what we just have said, the BCS state (4) still holds true at V 2 = 0. 
(However, for V 2 = 0, the critical temperature T c = 0, see Sect. 3). Another way to 
see this stems from the usual Thouless criterium for the onset of superconductivity 
which we here consider using the Bethe Salpeter equation for the 2-particle Green 
function (in a symbolic writing) [4] 

G 2 = G° + G°V 2 G 2 , (20) 

where G° = (1 - 2n p )/(u - 2(e p - //)). For attractive V 2 ^ eq.(20) yields the 
usual Cooper pole singularity. However, in systems with FC, eq.(20) shows a pole 
at uo = for the condensate state with e p — /j, even if V 2 = 0. 

A consequence of these considerations concerns the entropy S of systems with 
FC. Since the ground state wave function (4) does not change its BCS structure as 
a function of V 2 down to V 2 — and since \BCS > is a pure state, we conclude 
that the entropy S of the system is zero even at V 2 = 0. Thus, the entropy paradox 
inferred from a nonzero value of S(T = 0), when calculated on the basis of the 
ordinary Fermi gas formula 




with the distribution (12), is removed. Instead we should use the BCS expression 
for S (see below) which again shows that S = at T = 0. 

The superfluid properties following from the formulas obtained, turn out to be 
quite unusual. First, in the general case of finite range forces when the strength 
of the pairing term V 2 << V\, we seemingly deal with the weak coupling limit of 
BCS theory. However, in this first shown in [2] and also seen from (16), 

the gap in the spectrum of the single particle excitations has no BCS exponential 
smallness: A BC s ~ exp(— e° F /V 2 ), since on the contrary, it turns out to be linear 
in the particle-particle (p-p) effective coupling constant V 2 . Furthermore, outside 
the interval pi < p < pf, the gap A p vanishes, the boundaries pi and pf practically 
depending only on the value of V\. Thus, the gap A p exists only in the region 
occupied by the Fermion condensate. These properties and the fact that k p ^ 
even for V 2 = (and thus A p = 0) show that the strong BCS pairing is triggered 
by V\ (in conjunction with the quasiparticle plateau). Pairing with such unusual 
properties is a "shadow" of the Fermion condensation. 
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A quite unexpected point is that a nontrivial solution A p 7^ exists in our 
frame of singlet pairing even if V 2 becomes negative, i.e. repulsive, contradicting 
completely the Cooper pairing pattern. This can be seen in the following way. 
First of all, nothing dramatic happens if the sign of V 2 in A p (10) and k p (1A, 14') is 
turned around. Of course, the form V 2 n p in the energy (6) indicates a loss of energy 
for finite k p rather than a gain as this was the case for V 2 > 0, i.e. attraction. 
However, the other option, which is no pairing, means unambiguously that n p is 
back to the Fermi step np(p) = 9{p — Pf)- Because of the square of the momentum 
distribution n p this in turn entails an energy loss in the term Vin p even greater 
than if we had kept k p finite (V 2 repulsive!). We convinced ourselves from numerical 
studies that the inclusion of the kinetic energy does not turn over the situation of 
the energy balance, at least for the case \V 2 \ << V\ considered here. That a gain 
in energy prevails even if V 2 is repulsive can also be seen from eq.(18): the sign 
of the correction to the chemical potential — (14 + V 2 ) 2 /4&e F is independent of the 
sign of V 2) as long as \V 2 \ « V 1 ( otherwise higher order terms must be included). 
Another way to characterize the situation is that FC is analogous to Bogoliubov 
pairing in a nonideal Bose gas. Indeed, according to (10), V 2 is nothing but the 
factor of proportionality in the relation between the gap A p and the anomal density 
K p . In contrast to BCS theory, k p is practically independent of the p-p interaction 
V 2 « V\ (see eq.(14),(14')). It is "prepared" by the interaction V\ entailing a 
non- Fermi-step distribution n p (see Fig. 2). 



3 Finite temperature 

The extension of our model to finite temperature is, in principle, straightforward. 
It follows along the lines of ordinary BCS theory [4]. First, the quasiparticle distri- 
bution n p and the anomal density k p are modified as follows 

n p = v 2 p (l - f p ) + u 2 p f p , k p = u p v p (l - 2f p ), (22) 

where f p =< a p a p > is the distribution function for the quasiparticle excitations 
(see below). They are described as a gas and the entropy S is given by the integral 

S = ~2 J 7^[/ P ln/p + (l-/ P )ln(l-/ p )]. (23) 
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We have now two different equations of the minimum: = and Ip- = with 
VL = E — fiN — TS. The last one is written as 

&(1 - 2v 2 p ) + 2u p v p A p - T In = 0. (24) 

where (e p — fx) — £ p as before. When deriving this equation we again have used the 
definitions e p = j^-, and A p = — j^-- The first equation of the minimum keeps the 
same form as at T = (eq.(8)). Therefore introducing the quasiparticle energy E p 
by a finite temperature generalization of (16) 

E P (T) = ^ 2 P (T) + Al(n (25) 

we can, as before, solve eq.(24) yielding the ordinary relations for the u- and v- 
factors: 

Then inserting these relations together with (25) into (24) we are led to a seemingly 
ordinary Fermi-Dirac distribution 

fp = ■ (27) 

e t + 1 

As a matter of fact, this relation represents an equation for f p since E p implicitly 
depends on f p . 



a. The case T < T c . 

Let us turn to the solution of the equation for the gap A P (T) following from its 
definition and the relation (22) for k p (T) 

A P (T) = V 2 u p (T)v p (T)(l - 2f p (T)) (28) 

Upon substituting into this equation relations (25-27) we get 

We shall see that a nontrivial solution of this equation written as 

2£ p (T) = V 2 tanh^|^, Pl (T) < p < p 2 (T) (30) 



only exists in an interval [pi,j?2], narrower than the interval [pi,Pf] (13) of FC. Being 
independent of p, E P (T) drops with T and vanishes at T — T c — 14/4 [3]. This is 
also seen in Fig. 3 where a graphical solution of (30) is given. 

From (30) one finds tanh |=p = 2E p /V2 which allows us to rewrite (27) as 



f P (T) = \-^, Pl <P<p 2 . (31) 



With this relation, eq.(22) for n p takes the same form as the first of eqs.(12) with the 
only replacement of fi(T = 0) by n{T). Taking into account that the second relation 
between these quantities following from the definition e p = 4^- stays unchanged (see 
(9)) we have for n p and e p practically the same form as at T = 0: 

f x(T) + V 2 /2-p 2 /2m 

'W ) - v + v ' 

_ ^ = MiM^ P1<p< P2 . (32) 

V\ -+- V2 

Outside this interval the gap A p = and the quasiparticle distribution becomes 

1 . . 

n P = ep( T)- M (T) > P<Pi or p > p2, (33) 

e t +1 

as can be seen from the first of eqs.(22). Since relation (9), linking the quasiparticle 
energy e p to the distribution n p , also holds at T ^ 0, this means that eq. (33) is 
a quite intricate equation rather than an explicit formula for the calculation of n p . 
We can rewrite relation (33) between these quantities in the form 

e„-/i = Tln^^ (34) 
n p 

and substitute it into (9). Then one obtains 

1 — 7% 

p 2 /2m — fx + V\n p = e p — /i = Tin -, p < pi, or p > p 2 , (35) 

Tip 

This is a closed equation for the evaluation of n p . Finding n p we can insert it into 
eq.(34) to calculate e p . This procedure will be considered in more detail for the case 
T >T C later on. 
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Now we turn to the calculation of the boundaries pi(T), P2CO, determining 
the volume occupied by the superfluid phase. They are found from matching the 
distribution n p given by the formulas (32) and (33). At these boundaries, eqs.(32) 
and (33), we can rewrite the matching equation in the form 



or equivalently 



*^- 2 (36) 

2£(^ = \/ 2 tanh^ (37) 

where pk = P\ or p 2 . It is the same equation as (30) but for A p = 0, as it must be. 
It also has a nontrivial solution only if T < T c — V 2 /A. The distance between the 
points Pi(T) and p 2 {T) determining the volume occupied by the superfluid phase 
diminishes with T and they go to meet each other at T = T c at the midpoint pu of 
the interval [pi,Pf] where n PM = 1/2. 

Upon expanding the r.h.s. of (37) which vanishes at T = T c one finds for the 
support of the superfluid phase close to T c : 



The same procedure applied to eq.(30) yields 



£ p = 27^3(1-^), Pl <p<p 2 . (39) 

Thus, the gap in the spectrum of the single particle excitations has a quite unusual 
behavior. It is limited by the boundaries p\ and p 2 and the support shrinks as T 
approaches T c . Since the absolute value of the gap also diminishes similar to the 
usual BCS gap, it is like "an ice cube which is melting". This behavior is shown in 
Fig. 4. It can easily be obtained from A p = ^ E^{T) — Q(T) using the solution of 
eq.(30) shown in Fig.3. It is in strong contrast with the T-dependence of the usual 
BCS gap which "melts" only in "height" but whose p-extension, given by the Debey 
temperature, is independent of T. The shrinkage of the volume of the superfluid 
phase near T c significantly influences the density of states p(e) and may exhibit itself 
in tunneling phenomena. 
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b. The case T > T c . Let us now turn to the normal state of the system, 
i.e. to T > T c . Since then A p = k p = 0, eqs. (33,34 ) and (35) hold in the 
whole domain. It is straightforward to solve eq.(35) with (33) for e p . The result 
is shown in Fig. 5 for T << e° F (T c can always be chosen arbitrarily small in the 
case V 2 « V± « e° F considered here). We see that in comparison with the zero 
temperature case of Fig.l the plateau is just a little tilted and rounded off at the end 
points. Inserting e p (T) of Fig. 5 (see also eq. (41) below) into (33) we realise that 
even for T > T c the distribution function has qualitatively the same behavior as the 
one seen in Fig. 2 with only corners rounded by the temperature effects. The spread 
of n p for T > T c is therefore still governed by features of the Fermion condensation 
and not by temperature. Thus at T c < T < V\ we are faced with a quite interesting 
situation. Superfluid current is absent while phenomena which occur due to the 
spread of the momentum distribution n p , usually associated with superconductivity 
(e.g. deviations from the Korringa law in the spin-lattice relaxation rate [5]), persist. 

To have more quantitative and analytical insight into what is going on let us 
first find the low temperature expansion for n p in the limit V 2 — > 0. For small T the 
r.h.s. of eq.(35) is small. The lowest order of n p is obtained in omitting the r.h.s. 
completely. The result, of course, coincides with (17). The next correction is 
found by the substitution of (17) into the r.h.s. of (35) that yields 

(1) _ /x-pV2m T K-y+Wm 

This correction contains a small parameter T jV\. The iterations can be continued 
giving rise to a low temperature expansion of the quasiparticle distribution n p in the 
case A p = 0. There is a one to one correspondence between iterating the solutions 
for the distribution n p (T) and the spectrum £ P (T). To linear order the latter is given 
by (see eq.(34)) 

1 - 4 0) rp, Vi-^ + p 2 /2m 2 
= Tln = Tln n-jfifim + 0iT ] - (41) 

In the limit e p — > fj,, at T finite, we can further simplify (41), since then the argument 
of the logarithm must be close to one. This leads to 



£ P (T) ~ 2T 



K 



+ 0(T 2 ). (42) 



V 1 

This result can also be obtained more directly upon inserting (33) into (35) and 
taking the limit £ p — > at finite T << V\. However (42) is valid only in the 
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immediate vicinity of e p to /i whereas (41) holds for a longer range of p- values. 
From (41) and (42) we see that the quasiparticle spectrum e p (T) in the region of the 
Fermion condensate depends linearly on T. This feature is very different from the 
ordinary Fermi liquid with its T 2 dependence of e p (T) stemming from the Landau 
formula for the variation of the quasiparticle energy 



5e p (T) = 2j^f( P , Pl )5n Pl , 



where / is the Landau scattering amplitude. A linear temperature law may seem 
quite surprising, since in the case considered the specific heat Cy = OE/dT = 
TdS/dT = -TZ p (e p - fi)dn p /dT is a linear function of T [3] and hence E(T) - 
E(0) ~ T 2 as in ordinary Fermi liquids (this is easily understandable, since for FC 
we have £ P /T = const (see eq.(42)) and thus dn p /dT = 0. For this reason the 
condensate does not contribute to Cy). 

Let us therefore pin down the point where the conventional argument leading to 
the Independence of e p (T) at low T fails. To this purpose, we differentiate both 
sides of the previous equality with respect to T and obtain 

de p (T) f d 3 Pl f pi 

~W~ i (2^p /o(p ' Pl) r? n f l(1 ~ n ^~ (43) 

Here f is the s-wave part of / and use has been made of (33) for dn p /dT = 
(s p — p)n p {l— n p )/T 2 neglecting, as usual for normal Fermal liquids, the insignificant 
T-dependence of the quasiparticle energies on T. One can then introduce x p = £, P /T 
as a new integration variable and extend the integration limit at low T with good 
accuracy to ±oo. Since xn(x)(l — n(x)) is an odd function of £ p we have to expand 
the remainder of the integrand, i.e. , li /e p ~/(de p /dp) to linear order in £ p and hence: 
de p (T)/dT ~ T or e p ~ T 2 . In a system with FC this demonstration, however, 
does not hold, since (de p /dp)~ l diverges exactly at the Fermi surface. Therefore all 
arguments in favor of the T 2 behaviour of the spectrum e p {T) fail and formula (41) 
survives. 

A straightforward comparison demonstrates a remarkable correspondence be- 
tween the present results and those of the Landau theory: most relations of ordinary 
Fermi liquid theory also hold for systems with FC, we simply have to remember that 
the effective mass diverges as T _1 : 

m = ;;Ur)„ = 27V, -,?) = m 3r (44) 



12 



We will use this relation in Sect. 4. (For quantities which contain the derivative 
dn p /dT like the specific heat Cy ~ pfiti*T, the insertion of (44) is not valid, since 
otherwise Cy would be a constant in T in contrast to what we have seen above). 

The memory of FC for temperatures beyond T c can be seen in a particularly 
transparent way from the density of states 

Ke) = 2f-^S(e-e p {T)). (45) 

This integral can be calculated analytically since de p /dp can be obtained by differ- 
entiating eq.(35). The result is 

p(e) = ^(2m) 3 / 2 [ £ -y 1 n( £ - /1 )] 1 / 2 [l + ^( £ -^- £ )] 

- n °f f 2mT l n(e - p)n(p- e), n ~ e, (46) 

where n(£) = (e^ T + 1) _1 and N F = ppm/n 2 . For typical situations T c < T < V\ 
we show p{e) in Fig.5. 

We see that the condensate contribution is spread near the Fermi surface over 
the interval 5e ~ T and it is enhanced by the factor (pj — pf)/2mT ~ V±/T > 1 
if T ~ T c . Its magnitude drops as T _1 . It becomes of the order of the regular 
contribution N F to the density of states at a characteristic temperature, denoted 
further by Tf which is 

Pf—p iQ Q c 

~ 2^T~ ~ (47) 

where f2 c stands for the condensate volume and Qp, for the volume of the Fermi 
sphere. Above Tf the effects of FC become insignificant. With (47) we can write 
the enhancement factor p{e = p)/N F of the density of states near the Fermi surface 
in terms of Tf as 

p(e = n)T f 



N° F T 



(48) 



It is worth noting that in such a form the estimates (47,48) are still valid for an 
arbitrary form of the effective particle-hole interaction as long as it gives rise to FC. 
We also realise that there is not a definite T value at which the effect of FC stops, 
its influence diminishes gradually. 

As a further peculiar point characterizing FC we want to consider the entropy 
at T > T c . Using (23) in the low temperature limit we obtain 



13 




d 3 p 



-n p (T = 0) In n p (T = 0) ~ V x In 




(49) 



(2k) 



Thus, the entropy of the system with the Fermion condensate at T > T c drastically 
exceeds ordinary superconductor values Sbcs(T c ) ~ T c ~ V^. This is a specially 
peculiar "visit card" of FC in homogeneous systems. 

Summing up this section we may say that systems with the Fermion condensate 
exhibit a quite unusual non Fermi liquid temperarture behavior. The main charac- 
teristic features are that the g £IP CIS db function of T not only shrinks in magnitude but 
also in phase space volume and that the level density enhancement and flattening 
of the quasiparticle spectra, i.e. the typical Fermion condensate features, persist up 
to temperatures T ~ Tf far greater than T c . Furthermore, the entropy for T > T c 
is drastically enhanced with respect to ordinary Fermi liquid values. 

Let us just add a word of the behavior of the various quantities as a function of 
T for T < T c . It is readily realized that the presence of a finite gap barely alters 
the characteristic features found for T > T c indicating once again that in this kind 
of physics FC is the driving mechanism and not pairing. The latter is a host, itself 
strongly influenced by FC. 



All the results of this article have been obtained on the assumtion that the approach 
based on the quasiparticle picture is valid implying that the width 7 of the relevant 
single particle states (in our case, the condensate) does not exceed the quasiparticle 
energy e p . This assumption is fulfilled for superfluid systems with the Fermion 
condensate because of the presence of the gap in the spectrum of the single particle 
excitations. Thus, aXT <T C the quasiparticle picture is valid. 

We shall see that the width 7 remains small even if T exceeds T c . This im- 
plies that the quasiparticle picture survives turning to the normal state of these 
systems. Considering this problem we can treat it, as it was said before, within the 
ordinary Fermi liquid approach by simply taking into account the T dependence of 
the effective mass m*(T) ~ T _1 (see eq.(44)). Much work has been done within 
the ordinary Fermi liquid approach to simplify the calculations and we shall make 
use of well known results. First of all, the width 7 is inversely proportional to the 



4 Validity of the quasiparticle picture 
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collision time for which the following formula can be derived [6-8] 



1 (m*) 3 T 2 W(9,<p) 

7 = < >' (50) 



cos 



2 



where W(9, 0) is the transition probability depending on the angle 9 between the 
vectors pi and P2 of the incoming particles and the angle between the planes 
defined by the vectors Pi,P2 and p^p^- The brackets < ... > denote angular 
averaging. 

It is useful to introduce a dimensionless quantity w defined as follows: < W/ cos | >= 
2txw/N f where N F = p F m* /n 2 is the density of states at the Fermi surface. With 
this result the formula for the collision time To becomes 



1 nm*T 2 

' ■■— — w > ( 51 ) 



where [6-9] 
w — — [ d(f) f d cos 6 

2,71 Jo Jl 



To 4 p F 



1 \A s {9,(t)) + A a {9,(t))\ 2 1 \A s (9,cf)) - A a (9,(f>)\ : 



4 cosf 2 cosf 



, (52) 



and A s and A a are components of the scattering amplitude. The triplet A T and 
singlet As scattering amplitudes are related to A s and A a as follows: At = A s + A a 
and As = A s —3A a . Let us note that in a Fermi liquid, without spontaneous breaking 
of symmetry, the knowledge of the component A s (9,4>) is enough to reconstruct on 
the basis of symmetry relations the other component A a (9,(f)). These relations are 
written for the triplet and singlet scattering amplitudes separately 

4r(Pi» P2> P3, P4) = -^t(P2, Pi, P3, P4), 

A S (Pl, P2, P3, P4) = ^s(P2, Pi, P3, P4)- (53) 

Substituting here A s and A a leads to 

A a (cos#,g 2 ) = -^A s (cos9,q 2 ) - ^A s (cos9,q' 2 ). (54) 

We have used that cos 9 remains unaltered under the interchange of the initial mo- 
menta while the transferred momentum q is replaced by q' = p 3 — p 2 . One has [8]: 

q 2 = (p 3 _ pi )2 = 4p | sin 2 I gin 2 I and q /2 = 2 p 2 F (l- COS 9) - Q 2 . 
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The variables cos 6 and q are particularly convenient if A s is linked to the Landau 
scattering amplitude / at an arbitrary value of the momentum transferred [6,9]. The 
corresponding equation has the form 

A'( P i, pa, q) = /( Pl , P2 , q) + 2 / ^/(Pr, p, q) A'(p, p 2 , q). 

(55) 

This equation is, as usual, solved by expanding both sides in Legendre polyno- 
mials. The familiar result for each component is: 

A s (a) - fl{q)NF (56) 

A(q) "l + iV^(q)L(q)/(2/ + l)- (56) 

Here the dimensionless Lindhardt function L(q) = —Xo/Np is introduced where 
Xo(q) is the linear response function of noninteracting quasiparticles 

X,(q) = 2/|^ -WW , (57) 

7 (27r) 3 £ p+q/2 - Ep-a/2 ~ V 

where the energy transfer u = e p > — e pi . This function exists in closed form [6] and 
varies smoothly with the parameters q/pF and oo/qvp. 

Since nothing is known about the Landau parameters // with I > 2 we keep only 
the 1=0 and 1=1 harmonics of f(6) neglecting all others with I > 2 as is usual in 
Fermi liquid theory [9]. In addition, we shall consider fo(q) and fi(q) as smooth 
functions of q. Hence, for T — > we shall neglect the "1" in the denominator of 
(56), since the density of states diverges for T — > 0. Then the scattering amplitude 
A can fully be constructed. As a result, one obtains: 

JN l + 3cos6> 

tarn i\ 1 + 3 COS ^ 1 + 3 COS 6 

A(eA) = -^L^T' 2 ^(V (58) 

Thus, we infer that in the strong coupling limit the scattering amplitude A(9, <ft) 
becomes a universal function independent of any parameter characterizing the in- 
teraction between particles. 

At this point one may wonder whether screening corrections to the Landau scat- 
tering amplitude / in (55) do not play a very important role as well (this point 
had also been stressed in [3]). In fact, this is not the case. A convenient frame to 
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discuss this point are the selfconsistent equations for / derived by Babu and Brown 
[10]. In that work / is split into a direct and an induced part. The direct part is 
a Bruckner G-matrix and the induced part f md corresponds to an exchange of an 
(RPA) phonon between the particle and the hole (schematically: / = G + f md ). For 
our purpose it is sufficient to restrict to the s-wave part f of /. Then the induced 
part is given by (eq.(4.23) and eq.(5.20) of [10]): 

find _ 9 f 2 ^ F L(q) 

As before we can neglect the "1" in the denominator of (59) and obtain f md = — 2/o 
which is a finite result. Therefore screening corrections do not invalidate our above 
conclusion. 

As mentioned already, the g-dependence of the Lindhardt function is insignificant 
and, in first approximation, it can be ignored by taking L(q) = L(0) = 1. Then w 
is evaluated analytically and after a simple integration one obtains 

w = —. (60) 
5 

With this result the width 7 is calculated in closed form: 

1_ l6nm*T 2 
r 5 p 2 F 

Inserting here m*/m = T//4T from (44) we finally find 



7C0 = - = — -ZS-- (61) 



7 (T) = ^T. (62) 

Thus, 7(T) is a smooth function of T, the width of the condensate states turns out 
to be rather small, and no inconsistency emerges in our calculations. We therefore 
can infer that the quasiparticle picture turns out to be valid in all the temperature 
interval < T < Tf of our interest. 

It should be mentioned that, had we treated 7 in the Born approximation, the 
quasiparticle lifetime would go to zero with T as in [3]. The reason for this difference 
stems from our interaction which takes into account screening effects. In different 
words, within perturbation theory the product of the effective interaction and the 
density of states goes to infinity. However, if screening effects are properly taken 
into account this product stays of the order of one in these systems. 
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It is useful to discuss also the e- dependence of the imaginary part ImE(p.e) of 
the mass operator S(p, e) in systems with FC. Usually, to obtain the ^-dependence 
of ImE, if its T-dependence is already known, it is enough to replace T by £ = 
e — /i. This rule implies that in the system with FC, the imaginary part of the mass 
operator is proportional to £ in contrast to the £ 2 - dependence inherent in normal 
Fermi systems. A straightforward calculation does confirm the linearity of ImE as 
a function of £ [2]. The difference stems from the fact that, in systems with FC 
one of the three quasiparticles of the decay channel can belong to the condensate. 
Since the condensate energy is fixed, there is no additonal integration related to this 
particle and the phase volume has a value as if only two-quasiparticles were present 
in the decay. This leads to the linear behavior of ImE(p, e) in e — pi. 

5 Discussion and Conclusion 

In this article, we have again carefully analysed the implications and the physics 
of the Fermion condensation phenomenon. This phase transition was postulated to 
occur in strongly correlated Fermi systems several years back [2]. Our investigations 
were carried out in a slightly generalized version of the model of infinite range forces 
[3] which enabled us to analytically demonstrate the non Fermi liquid behavior of 
systems with FC. In our case we decoupled the link V\ = Vi considered in [3] between 
the parameters of the particle-hole and particle-particle channels and analysed the 
case V 2 « V\ « e° F . The first thing which our study revealed is that the BCS 
ground state persists even if pairing forces (V^) are absent. The smearing of the 
quasiparticle momentum distribution n p is therefore due to the particle-hole (Vi) 
rather than the particle-particle (V2) interaction. This implies that the system with 
FC is simultaneously in the paired state with phase locking independent of whether 
there is a pairing force or not. Since the BCS ground state is a pure state there 
is no residual entropy problem at T = as this was evoked earlier [3]. Even more 
dramatic, we found that, due to an unusual balance, the BCS solution persists even 
for repulsive Vi as long as | V2I << V\. In these limits the BCS state therefore 
exists independent of the sign of Vi and the entropy always is zero as T — > 0. At 
finite T there exists a critical temperature T c = V 2 /A where the gap disappears. 
This disappearance occurs, however, in a non-Fermi-liquid fashion: A p not only 
shrinks in magnitude but also in momentum extension (phase space) which may be 
a measurable effect. Once the gap has disappeared the system still shows Fermion 
condensation behavior up to temperatures Tf ~ V\ > T c : the typical flat plateau in 
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the single particle spectrum e p is only slightly tilted (~ T) and rounded off at the 
edges and the fall off width of n p is of the order of V\ and not of the order of T < V±. 
This is another example of non-Fermi-liquid behavior. The effect nicely shows up 
in a strong enhancement of the level density at the Fermi level for T > T c < Tf. 
One can also say that the standard Fermi liquid expressions remain valid with one 
replacement: for T > T c the effective mass behaves itself as m* ~ T _1 . 

An important issue of FC is the validity of the quasiparticle picture. This was put 
into question in [3]. However, there particle collisions have been treated in the Born 
approximation whereas the diverging density of states demands a nonperturbative 
treatment. Including polarization effects yields a quasiparticle width 7 of the order 
of T and the quasiparticle concept is preserved. 

Of course, the phenomena analyzed here can occur not only in systems with 
infinite range forces. Solvable Hartree-Fock-like models suggested in [2,3] do con- 
firm the conclusion that the case of finite range forces merely changes quantitative 
details without changing qualitative features of the Fermion condensation such as 
the plateau in e p etc. In particular, the shrinkage of the volume occupied by the 
superfluid phase also survives with only one modification: the gap A p acquires an 
exponentially small tail instead of a sharp cut off. No change in the conclusions also 
occurs if we go beyond the Hartree-Fock approximation. For example, we can add 
to the functional (7) a term of non HF nature of the form SF ~ E PiPl n^n^. With 
this term the problem again can be solved nearly analytically and we can verify that 
all qualitative results, obtained above, hold in this model, too. The same is true 
for the relation (44) between the characteristic temperature Tf and the size of the 
region occupied by the Fermion condensate. Further, when calculating the width 
7, we did not use the HF approach at all so that the quasiparticle picture persists 
beyond the HF approximation. 

A separate question concerns the possibility of the existence of Fermion conden- 
sation in real physical systems. First of all, it is related to the evaluation of the 
single particle spectra in such systems. A theory of these spectra in normal Fermi 
liquid based on a functional approach has been constructed in [2] . It is in agreement 
with the RPA results obtained by Gell-Mann-Brueckner-Galitskii [11,12]. Results of 
calculation of the single particle excitation spectra carried out for electron systems 
within the jellium model demonstrate that the point in the dispersion e p , where 
the derivative de p /dp = 0, a precursor of the Fermion condensation, does emerge in 
these systems at the densities at which the effective interaction between electrons 
becomes strong enough [13]. In this case one can also calculate the width 7 for the 
condensate particles and again find j(T) ~ T validating the quasiparticle picture. 
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An issue of considerable interest is the finding [14] within the local density 
approximation (LDA) of bifurcated saddle points in double-plane materials like 
YBa2Cu3C>7, YBa2Cu4C>8, and E^S^CaC^Os+s. It seems that, for a certain range 
of doping, these saddle points are pinned to the Fermi level and that they lie ex- 
actly where experimentally an anomalous flattening of the single particle dispersion 
in two-dimensional strongly correlated electron systems is observed [15]. With the 
Fermi level in the region of the bifurcated saddle point again the Maxwell construc- 
tion seems applicable [16], leading to straight segments of Fermi lines instead of the 
two-dimensional plateaus discussed in this work. It should, however, be noted that 
these flat Fermi lines also lead to a diverging density of states pie = /i) [17]. We 
intend to investigate this problem in a future publication. 

To summarize, we have analysed the rather unusual properties of Fermi systems 
which undergo a phase transition to Fermion condensation. Strong non Fermi liquid 
behavior has been demonstrated for those systems whereas the the quasiparticle 
picture seems to survive: for T < T c due to the gap in the single particle spectrum, 
for T > T c , due to a screening suppression of the matrix elements responsible for 
the decay process. 
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Figure Captions 



Fig.l. Schematic view of a quasiparticle dispersion with a downswing at the 
Fermi energy and the corresponding Maxwell construction. 

Fig. 2. Normal (n p ) and anomalous (k p ) occupation numbers in the model with 
the Fermion condensate. 

Fig. 3. The quasiparticle energy E p (eq.(30)) as a function of temperature (in 
units of Vi). 

Fig. 4. The gap function A p = ^ — ^ as a function of £ p = e v — /x for various 
temperatures. 

Fig. 5. The quasiparticle dispersion around the chemical potential for various 
temperatures T > T c . 

Fig. 6. The level density p(e) corresponding to the same situation as in Fig. 5. 
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